February 6, 2008 



(N 



O 

o 



Quantum dissipation theory of slow magnetic relaxation mediated by domain-wall 
motion in one-dimensional chain compound [Mn(hfac)2BNOH]- 



A.S. Ovchinnikov* , I.G. Bostrem, V.E. Sinitsyn, A.S. Boyarchenkov, N.V. Baranov 
Department of Physics, Ural State University, 620083, Ekaterinburg, Russia 



O 
O 

' K. Inoue 

, I Department of Chemistry, Faculty of Science, Hiroshima University, Higashi-Hiroshima, 739-8526, Japan 

, ^ ■ Based on a quantum dissipation theory of open systems, we present a theoretical study of slow 

' dynamics of magnetization for the ordered state of the new molecule-based magnetic complex 

[Mn(hfac)2BNOH] composed from antiferromagnetically coupled ferrimagnetic (5/2, 1) spin chains. 
Experimental investigations of the magnetization process in pulsed fields have shown that this com- 
pound exhibits a metamagnetic AF-FI transition at a critical field in the order of the interchain 
coupling. A strong frequency dependence for the ac-susceptibility has been revealed in the vicinity 
' ^ ■ of the AF-FI transition and was associated with an AF-FI interface kink motion. We model these 

ly-j ' processes by a field-driven domain-wall motion along the field-unfavorable chains correlated with a 

dissipation effect due to a magnetic system-bath coupling. The calculated longitudinal magnetiza- 
tion has a two-step relaxation after the field is switched off and are found in good agreement with the 
experiment. The relaxation time determined from the imaginary part of the model ac-susceptibility 
agrees qualitatively with that found from the remanent magnetization data. 



I. INTRODUCTION 



, Magnetic resonance effects due to switching of spins by a time-dependent magnetic field and relaxation measure- 

Q • ments are widely used to study magnetic materials. They may give information about the main mechanism of the 

O ' magnetization change in samples. When a magnetic field varies, magnetic materials which exhibit a hysteresis present 
a characteristic time dependence of the magnetization due to the multiplicity of available metastable states. In many 

' relaxing magnetic systems the time dependence can be described by M = Mq — Slnt, where S is the magnetic vis- 

J> ■ cosity of the system. The viscosity is a consequence of thermal activation of irreversible domain processes such as the 

' domain wall motion and the nucleation of domains of the reverse magnetization. The logarithmic change with time 

: M ~ In i is predicted if there is a distribution of energy barriers or time-dependent activation energies present in the 

• material. A simple Debye relaxation M = Afoexp(— t/r) arises from a single-barrier activation mechanism^'^. The 

^ , barriers associated with the relaxation process are of two types. The intrinsic barriers arising, for example, from the 

' magnetic anisotropy contribute to the reversal of magnetization, whereas the barriers due to the pinning of domain 

[ walls are generally attributed to the defects in the materials. Both types of barriers are responsible for the pronounced 

"Y*^ ■ metastability (hysteresis phenomena) of magnetic systems. 

' Recent efforts in synthetic chemistry provide a number of low-dimensional magnetic systems that show the slow 

, relaxation of the magnetization, for example, this effect was found in one-dimensional (ID) anisotropic ferrimagnetic 

I ' chains named as single chain magnets (SCM)^. The slow magnetic relaxation in the paramagnetic phase has been 

' ^ , observed with ac susceptibility and SQUID magnetometry measurements in the real quasi-ll? ferrimagnetic compound 



[Co(hfac)2NITPhOMe] . An Arrhenius behavior with the activation energy A ~ 152 K, which is of order of the 



intrachain exchange interaction between alternating Co^+ and organic NITPhOMe spins, has been observed for 10 
decades of relaxation time and found to be consistent with the Glauber model^. The relaxation was also studied by 
. 5^ ' monitoring the decay of the longitudinal magnetization, which was found to be exponential. The slow relaxation is 
governed by the uniaxial anisotropy seen by each spin on the chain and magnetic correlations between the spins. 
Recently, the ac susceptibility and magnetization in steady and pulsed fields have been measured for the new 
5^ , molecule-based magnetic complexes [Mn(hfac)2BN0R] (R=H,C1) with ferromagnetically (CI) or antiferromagnetically 
(H) ordered ferrimagnetic chains composed of S—1 (biradical) and S=5/2 (Mn^+)^. It has been found that the change 
in magnetization in these compounds under application of a magnetic field below the three-dimensional (3D) ordering 
temperature (Tn =5.5 K (H) and Tc =4.8 K (CI)) is a slow dynamical process which presumably originates from 
their strong one-dimensional character, i.e. because of the weakness of the interchain exchange (J) interaction in 
comparison with intrachain (J') one (J/ J' ~ 10^^). The thermally activated change of the remanent magnetization 
in [Mn(hfac)2BN0H] after switching off the field looks quite different from that in [Mn(hfac)2BN0ci]- Figure 1 
displays time dependence of the remanent magnetization at 1.5 K for these compounds for comparison. The large 
value of the relaxation times of both processes and their existence well below lead to the suggestion that these 
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relaxation processes are related to the development and motion of magnetic domains either with a wide distribution 
of energy barriers (R=C1) or a single energy barrier height (R=:H). 

Indeed, the change of the magnetization of [Mn(hfac)2BN0H] with time during and after application of a pulsed 
field is controlled by the direct AF-FI and inverse FI-AF transitions. For the applied pulse with a duration of 20 
ms and amplitude of 4 T, the direct AF-FI transition is completed within a time of less than 2 ms. The high rate 
of the direct transition is due to the high value of the magnetic field in comparison with the critical field of the 
AF-FI transition (0.03 T in steady fields and about 0.2 T in pulsed fields). The large relaxation time (^500 ms) after 
removal of the pulsed field implies very slow dynamics of the magnetization during the inverse transition from the 
field-induced FX state to the initial AF state in zero field. The slow dynamics of the magnetization in this material 
was suggested to be be controlled by the domain wall motion along the separate field-unfavorable chains- . The energy 
barrier which hampers the reversal of the magnetization is originated by the antiferromagnetic interchain coupling 
that is reminiscent of the molecular clusters or single cluster magnet (SCM) where the energy barrier is due to the 
magnetic anisotropy. 

The fast relaxation during the field change and the slow relaxation due to the transitions between the different 
metastable states involves the release of the magnetic energy which cause local heating. The process of energy losses 
in [Mn(hfac)2BN0H] is likely due to the motion of domain boundaries within the chains and is connected with the 
coercetivity losses and the radiation of phonons. Finally, the energy of the magnetic system is transformed into the 
phonon energy, that is, into the heat. Due to the magnetic system-bath coupling we may explain the experiments 
under the condition of isothermal relaxation, and, as a consequence, in an assumption that self-heating in the magnetic 
relaxation does not destroy 3D order. 

II. PRELIMINARIES 

The choice of a proper formalism of the quantum dissipation theory depends essentially on the real physical media 
and real process of the energy transfer. The interaction between a quantum system and its environment is the 
physical factor responsible for the relaxation process in the system. Such a relaxation process can be provided by the 
environment if it acts as a bath. This arises the fundamental problem of description of the relaxation dynamics for a 
system simultaneously interacting with a heat bath and a time-dependent driving field. 

The question now arises as to what is a nature of heat bath. We note that the spin system in the external dc bias 
field in ac-susceptibility measurements or in the pulsed field in magnetization measurements consists of two weakly 
interacting parts, i.e. field-favorable (Si) and field- unfavorable (^2) chains. Then, two distinct models for the bath 
are assumed. The first one identifies the bath with a crystal lattice (phonon bath). In this case, the energy levels 
of the bath are populated according to the conventional Gibbs distribution. When a rapidly varying strong pulsed 
field causes the AF-FI transition the spin temperature of the S'2 part becomes much higher than the lattice and the 
S'l part temperatures (saturation process). After the field is off the whole spin quantum system is driven to the 
thermal equilibrium characterized by Gibbs's density matrix owing to energy exchange between the S'l, S'2 systems 
and the heat bath. Following the fundamental Bogoliubov's procedure of the contraction of description with the 
accompanying hierarchy of relaxation timesi we may conclude that the weak interchain interactions are related to 
long-time relaxation process with the characteristic time equals to the characteristic time scale of the experiment 
1 — 10^ ms. It is assumed that another interactions, including usual spin-phonon ones, are strong enough and have 
associated correlation effects with shorter relaxation times, i.e. much smaller than the duration of the experiment. We 
should eliminate them as an irrelevant information to characterize the macrostate. Apparently, a possible approach 
to this irreversible process may be performed within the nonequilibrium statistical operator (NSO) method largely 
elaborated by Zubarev and co-workers^, which is a large generalization of the Gibbs's theory. The results of the 
NSO analyses will be given by us elsewhere. 

The second model of the bath is based on the suggestion that the relaxation process is related with the development 
and motion of magnetic domains in the field-unfavorable chains S'2 in a driving field. The recent study of the behavior 
of domain walls in Ising ferromagnetic chains yields the following pictureiS. In zero field the interface between two 
domains of oppositely oriented spins, a kink, moves left or right with the probability p = 1/2, which may be interpreted 
as a random walk. Two such kinks can meet via diffusive motion. Once there is only one spin left between two kinks 
they annihilate with probability p ~ 1 in the next time step and the two domains merge. Switching on the external 
driving field causes domains with spins parallel to the field to start growing. The kinks at the end of such a cluster 
move outwards one step during each time step where the field remains favorable. If the field switches into unfavorable 
direction, the cluster shrinks again (breathing behavior). In this process the small domains and the kinks associated 
with them will be eliminated rapidly. In the Ising model a domain wall width is simply the lattice constant. On the 
contrary, the domain walls in the Heisenberg chains are much wider due to the strong exchange interactions involved. 
The closest analog of such wall is the soliton in a magnetic chainii. Therefore another model, appropriate for this 
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domain wall relaxation dynamics, may be applied for the bath. In this model a quantum system (the field-unfavorable 
chain with a kink) interacts weakly with the environment, i.e. with the nearest surrounding field-favorable chains 
and the lattice (Fig. 2). The bath is again considered as a quantum statistical system being at equilibrium. The 
domain walls oscillate around their equilibrium position under a varying ac field (the process far from a saturation) 
or move reversibly under a time-dependent pulsed field (the saturation process). The bath provides an existence of 
random fields (Langevin forces) created by the environment which interacts with variables of the quantum system 
(Langevin dissipative modes). Due to the system-bath interactions the domain wall dynamics becomes irreversible 
and the system relaxes from its initial nonequilibrium state to the equilibrium one when the field is off. 

A quantum dissipative theory (QDT) with the system-bath interaction being treated rigorously at the second-order 
cumulant level for reduced dynamics has been recently constructed for open quantum systemai^. The theory belongs 
to a class of widely used quantum master equations, such as the Bloch-Redficld theory and a class of Fokker-Plank 
(FP) cquationsSiiS,, and is valid for arbitrary bath correlation functions and time-dependent external driving fields. 
The QDT-FP formulation constitutes a theoretical framework for the present study of dissipative processes in the 
molecule- based magnetic complexes [Mn(hfac)2BN0H]- 



III. QUANTUM DISSIPATION THEORY 

The key theoretical quantity in a quantum dissipation is the reduced density operator p(t) = Trspxit), i.e. the 
partial trace of the total system and bath composite pT{t) over all the bath degrees of freedom. For a system dynamical 
variable A, its expectation value 

{Ait))^Tr[ApTit)]^Tr[Ap{t)] 

can be evaluated with the substantially reduced system degrees of freedom. Quantum dissipation theory governs the 
evolution of the reduced density operator p(t), where the effects of bath are treated in a quantum statistical manner. 
The total composite Hamiltonian in the presence of classical external field can be written as 

HT = H{t)+HB-Y,QaFa (1) 

a 

Here H{t) is the deterministic Hamiltonian that governs the coherent motion of the reduced system density matrix 
and involves interaction with an arbitrary external classical field h{t). The system is embedded in a dissipative bath 
{Hb) and the last term in Eq.([T]) describes the system-bath couplings, in which {Qa} are Hermite operators of the 
primary system and can be called the generalized dissipative modes. The generalized Langevin forces 

Fa{t) = e'^«*Fae~'^«* 

are Hermite bath operators in the stochastic bath subspace assuming Gaussian statistics. Without loss of generality, 
their stochastic mean values are set to {Fait)) ^ = 0, where (. . ) g denote the ensemble average over the initially 
stationary bath density matrix pb(0). The effects of Langevin forces on the reduced primary system are therefore 
completely characterized by their correlation functions Cab{t) = {Fa{t)Fi,{0)) g. They satisfy the boundary conditions 
C'a6(±oo) = 0, and the detailed-balance and the symmetry relations C*^{t) = Cab{t — ifi) = Cba{~t), where (3 = 1/kT. 
This admits the Meier- Tannor parametrization Cab{t) in terms of exponential functions^ 



Cabit >0) = J2 



with the adjustable parameters v'^, These parameters are in general complex but, for simplicity, we take to 
be real and positive. 

The frequency-domain symmetry relation reads as C*{,(a;) — Cba(w), and the detailed-balance relation in terms of 
spectral functions is Cf,a(— w) = e^^Cabi^)- Using the generalized bath interaction spectral density function Jab{^^) = 
Cba{-oj) - Cab{ijj) obeying the symmetry relations Jab{ijj) = -Jba{-oj) = JLi^) have Cab{uj) = Jab{ui)/ {e'^'^ - l). 

We will use the reduced Liouville equation, i.e. the equation of motion for the reduced density matrix p(t), in the 
partial ordering prescriptionii 

m ^ -m)p{t) - R{t)p{t), (2) 

which is characterized by the local-time kernel R(t). The Liouvillian L is the commutator of the reduced system 
Hamiltonian H{t) in the presence of external classical field 

L{t)A=\H{t),A], (3) 
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and the superoperator R(t) can be formulated in terms of the system-bath interaction. In the standard approximation 
of the weak-couphng hmit in which the system-bath interaction is considered only up to second order the dissipation 
term is 

Rit)pit) EE J2 [Qa, Qa{i)p{t) - , (4) 

a 

where Qa{t) is the non-Hermitian relaxation operator in the Hilbert space 

t 

Qait)=Yl J dTCab{t-T)G{t,T)Qk. (5) 

The Liouville-space propagator G{t, r) associated with L{t) is defined via 

i^G{t,T) ^ L{t)G{t,T). 

It can be defined in terms of the Hilbert-space Green's function G{t, r) via the relation for an arbitrary operator A 

G{t, T)A = G{t, T)A&{t,T), 

where we treat A in the left-hand side as a vector in Liouville space. 

The reduced system Hamiltonian in the presence of external classical field can be written as 

H{t) = H, + H,fit), 

where is the time-independent, field- free Hamiltonian, whereas Hsf(t) is the interaction between the system and 
the external classical field h{t). We further define similarly to Eq.Q the Liouville superoperators Lg and Lsf{t) 
corresponding to the reference Hamiltonians. The identity 

t 

G{t,T)^Gs{t,T)-t J dr' G{t,T')L,f{T')Gs{T',T) (6) 

T 

obtained from the definitions 



G(t, r) =fexpl-ij [L, + L,/(r')] dr' j , 



where the symbol T implies a time ordering, and 

allows us to separate the dissipation effects in Eq.® into the field-free part and the correlated driving-dissipation 
part. This yields 

t t 

Qa{t)^Q:~iY, J dr J dr' Cab{t-T)G{t,T')L,f{T')Gs{r' -T)Qt, (7) 

-oo r 

where the field-free contribution is time-independent and given explicitly by 

t oo 

Qi = Yl I '^^ ^-^(* - ^)^-(* - ^Y.h^ Cab{r)e-''''^Qb - Cab{Ls)Qb- (8) 

- oo 

The equation ^ together with Eqs. (l4l7l8p provide a prescription for obtaining the reduced density operator up to 
the second order in system-bath interaction. The underlining assumption is that the system-bath coupling is not 
strong enough, which makes the second order cumulant expression reasonable. It is known that this approximation 
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applies well to most dissipative systems in quantum optics, and to transport in mesoscopic systemsis. In the strong 
system-bath regime a special technique is required, which goes beyond the second order approximationi2i2£. 

One of the traditional approaches to treat the problem is based on the associated quantum master equations. It 
focuses on the relation among various matrix elements of the density operator in the time-independent i/-eigenstate 
representation and is well suitable for the finite systems. For larger systems (spin chains, for example) the number of 
many-particle states increases dramatically and we cannot generally to solve all the microscopic equations. However, 
we can describe the system by macroscopic variables (domain wall position, total magnetization of a chain) which 
fluctuate in a stochastic way. The Fokker-Plank (FP) equation arises as an equation of motion for the distribution 
function of the fluctuating macroscopic variables^!. The Eqs. (|2l4l7l8p will serve as a starting formulations for deriving 
FP equations for observables of the reference system. 



IV. PROCESS FAR FROM A SATURATION: DYNAMIC SUSCEPTIBILITY 

Magnetic systems exhibiting relaxation phenomena can be characterized by the complex ac susceptibility, x(w) = 
x' — ix , where the dispersion x' and the absorption x are frequently dependent. Before moving on to the technical 
details of the calculation, we mention briefly some of the experimental results in the ac susceptibility measurements 
for the [Mn(hfac)2BN0H] compound that support the domain-wall motion picture. 

A strong frequency dependence both for x'i^) and x {^) has been revealed in the bias dc fields of 0.025-0.03 T. The 
magnetization process in the field range 0.02-0.05 T is accompanied by a remarkable hysteresis (of about 0.012 T), 
and, in addition, a small remanent magnetization was detected after removal of the field. These features are indicative 
of a magnetic phase transition of the first order, which occurs through a mixed phase state, from antiferromagnetic 
ordering of the chain magnetic moments to their parallel alignment in the field- induced state^^. In the region of 
the metamagnetic transition, where the AF and FI phase coexist, the amplitude of a maximum of both x' and x 
decreases significantly with increasing frequency, especially in the frequency region from 1 to 50 Hz. From the large 
values of x(^) for w ^ 1 Hz in the vicinity of the AF-FI transition, we may suggest that excitation of domain wall 
motion by a small oscillating field occurs more effectively at low frequencies. 

The complex magnetic ac susceptibility can be explained within approach, where the magnetization M is controlled 
by the field-induced sideways motion of domain walls. In this case, the contribution of one domain wall to the 
susceptibility is 

_ dM _ dM dx 
^ dh dx dh^ 

and taking approximately M oc a; as the magnetization increases due to a wall displacement along the x axis, one 
finds X dx/dh^. 

A periodic domain wall motion caused by the external ac-field h{t) ~ Hq coswt is modeled by a well studied system, 
a driven Brownian oscillator (DBO), with the Hamiltonian 

Hs = n (^a^a + , 

where the number of oscillator excitations a^a corresponds to an instant magnetization. The dissipation coupling 
mode Q 

Q ^ ^{q+a^ + a) 

interacts with a stochastic bath force. Here, q+ = are complex numbers, and a (a^) are the annihilation (creation) 
operators of the oscillator with the frequency determined by the interchain coupling in an applied bias field, i.e. in the 
vicinity of the AF-FI transition, and we hold only a single Langevin mode Q in study. The operator fi = ^ (^a^ + a) 

interacts with the ac-field and describes periodic domain wall movement caused by h(t). After introducing these 
definitions, Eq.([7]) can be transformed as follows 



oo 




(9) 
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where the system-bath couphng response 

A(t) = / dr' C{t') 

is given explicitly as 



A(r) = 



Q+^rn -C T -C r 



Cm + «^ 



(10) 



with the aid of Meier- Tannor parametrization. As usual, a summation is to be made for the repeated index to. In 
the following calculation we choose q± = I for simplicity, i.e. the dissipation is the same both for left and right 
domain- wall displacements. 

Substituting pU]) into Eq.®, followed by some minor algebra, we get 



Qit) = 



{C,jn COS cut -\- cu sin cut) . 



The time-local dissipation superoperator is 

i?(t)p(i) = [q, Q{t)p{t) - p{t)Q^ (t)] = Rsp{t) + im [Q. P{t)] 
where Rs is the field-free dissipation 



RsP{t) 



,Q'p{t)~p{t)[Q 



The effective local-field correction, acting on the system via Q is 



(Cm cos Lot + uj sin U)t) 



and determined by the imaginary part of the bath correlation function C{t). Hence, the final QDT formulation is 



The static superoperator Rs is 



Hs - phit) + Qm,p{t) - Rsp{t). 



Rsp{t)^ Q,Q'p-p{Qy 



(11) 



where the field-free time-independent dissipation coupling mode Q'' is 

Q' = C{Ls)Q = ^ (C {-n) a + C{n) flt) = ^ a + a)) , 
V 2 V 2 

where k+ — C (— il) and ~ C ( fl) are defined by the bath interaction spectrum C (il). Using the results of Sec. Ill 
we have C (fi) = J {n) n (n) and C {-n) = J (fJ) [n {n) + 1], where n (f7) = (exp (/3f7) - 1)"\ Then the dissipation 
superoperator Rg has the conventional formulation 



fc_ 1 k* 

Rsp{t) = —aa^p — — [k- + k*J) a' pa + -^paa) 



H — ^a^ap — - (fc^. -I- fc^) apa^ + -^pa^a. 



After some simple algebra, we obtain 

Rsp{t) = J (O) [n {ft) + 1] Q {ata, p} - apa^^ +J{n)n {Q) Q {aa\ p} - pa}j 
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Using the differential representation for the Bose superoperator, one can convert master equation (|lip into Fokker- 
Plank equations (see the Appendix for details) for the Wigner function / 



where z^z* are the complex variables, hit) ~ h{t) — ^{t) and 7 ~ J{il). To derive differential equations for Weyl 
the complex plane. Supposing that / — > at \z\^ —> 00, we obtain the system 



symbols (a)^,,, (a)Iy of the boson operators a and we multiply Eq. (jl2[) by z ot z*, respectively, and integrate over 



'w 



{^n + l){a)^+'-^, (13) 



Using the coordinate a; = (a + a^)^ / \/2 and the conjugated momentum p^; = — i <^a — a^)^ / ^2 we get the equation 
of the damped harmonic oscillator, 

^+7^+-o^--f^Mi) (15) 

with Wq = + 7^/4, which is the dynamic equation for the Bloch wall2^. 

The magnetic ac susceptibility was measured within the frequency range from 1 Hz up to 1 kHz (slow varying 
ac-field). In assumption that the characteristic time scale of the experiment of order 10^"^ 1 sec is much greater 
than the characteristic relaxation times of the bath (tj <C C)i we obtain 

i.e. the effective local-field correction ^(t) depends on the incident field. 

Now we use Eq.([T5]) to evaluate the range of relaxation time r ~ I/7. We note first that without an applied 
bias dc-field the oscillator frequency is determined by the interchain coupling, whereas in the bias fields of the 
experiment, approaching a critical value of the AF-FI transition, f2 is reduced to much smaller values, when a leading 
contribution to the ac susceptibility results from the motion of domain walls separating AF and FI phases. To reach 
a consistency with the data on a time evolution of the longitudinal magnetization in pulsed fields (see Sec. V), we 
suppose ri <C 7 and consider a small ac-field frequency w ^ ~ 7- Then the solution of Eq. psp has a relaxation 
character that yields the expressions for x'(^) ^-nd X (i-^) in usual Debye form 



X'(t^) = 



(l-cr)17 1 

~l 1 + a;2r2 ' 



X W = 2 -I , 2 2 ' (16 

where the relaxation time r = ^1 /oj^. The maximum of the imaginary part of the ac susceptibility x ('-^) is reached 
at Wmax = — loq / ^ ^ ^ . Accordiug to the available experimental data^ Wmax ~ 10 100 Hz (T = 3-^3.5 K) that 
yields t ^ 10 100 msec for the small frequencies ~ 1 Hz. This agrees qualitatively with the t values found from the 
relaxation of the remanent magnetization (Fig. 3). 

V. SATURATION PROCESS: STRONG MAGNETIC FIELD PULSES 

The key moment distinguishing this case from previous consideration is the value and time-dependence of the 
external driving filed. The ac-field being of order 10~*T is weak in the sense that the system remains near global 
equilibrium at all times. This is not the case for a strong field (~ 5 T) changing fast in comparison with relaxation 
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to global equilibrium. A long time scale of a driving field H{t) (~ 10 ms) prohibits the normal evolution towards a 
Boltzmann distribution of states due to dynamical non-Markovian effects. This feature is intrinsically built into the 
QDT, hence we can similarly construct a FP equation to evaluate the time evolution of a longitudinal magnetization. 
We introduce the operator fi = a^a interacting with an external pulsed field 



h{t)^ ho sin (^^t) m - 9{t - T)] 



where T is the time period of the external field and other definitions are identical to that used for the ac-field. This 
choice ensures a saturation of magnetization by pulsed field measurement performed in a half-pulse regime. 
Unlike the ac-case the commutator contained in the integrand of Eq. ([T]) now becomes an operator 

H'^') f„in(T'-r)„ _ „-iO(T'-T) t 



L,/(r')G,(T' - t)Q = -h{T') [fi, Gs{t' - t)Q] = (e'^^(-'--)a 
Further simplicity arises from that 



and 



G(t, r) a = exp i i J [n - h (r')] dr' 



G{t, t) at = exp i -i J [Q - h (r')] dr' I a^. 



Thus, we have 



G{t,T')L,f{T')Gs{r' -t)Q 



hjr') 
V2 



^in(t-T)-ihog(ty) ^ _ g-iO(t-r)+i/ios(t,T')^t 



(17) 



where / h{T")dT" = hog {t, t') . Substituting ^7}i into Eq.© we obtain 



Q{t) = Q^-—JdT A(t) h{t - t) 





where the system-bath coupling response becomes an operator 

oc 

A(t) = / dT'Cir') 



^iilr' -ihog(t,t-T) ^ _ ^-iflr' +ihog(t.t-T) 



(18) 



The explicit expression for A(t) can be easily carried out as 



A(r)- 



Cm + 



via Meier- Tannor parametrization. 

The convolution in Eq. (fT8|) is simplified as 



dr A(r) h{t - r) = < 



ho j dr A(t) sin [f (i - r)] , t>T 



t-T 



ho J dr A(r) sin [f (t - t)] , < t < T 



0, t<0 



The dissipative mode is then defined as follows 

Q{t) = (k+ a + fc_ a^^ 
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where 



and 



Cm 



Cm 4" 



(19) 



/ dr sin [f (i - r)] e^f^^-C'-^e^^'^of^*'*-^), t > T 



J dr sin [f (t - r)] e'"^-«'"^e-^''«9(*'*-^\ < i < T 



(20) 



0, i < 



The coefficients k+ = J (17) [n (fi) + 1] and — J (il) n (il) are determined as for the ac-case. The dissipation 
superoperator ^ reads as 



The Liouville equation 



takes the final form 



k 1 /- - \ k* 

R{t)p(t) — —^aa^p — — [k,^ + k*_j pa H — -^paa^ 



-—^a^ da — — ^fc+ + fc^^ apa^ H — ^pa^a. 



p{t)^-i[H,-fih{t),p{t)]~Rp[t) 



k 1 /- ~ \ k* 
p{t) — itl [p, a^a] = — —aa^p + — ^fc_ + fc* j a^ pa ■^paa'' 



—a'ap + - + apa' ^pa'a, 



(21) 



where il — fl — h{t). 

Then we convert Eq. pTjl into the equivalent Fokker-Plank equation using the Wigner functions for the density 
matrix p and Bose operators 



dn 
'dt 



^ (fc+ + fc;-fc_ -ki 

v2 



1 



1 



n + - I H \ k+ + k'^ + k- 



2V2 



fc* , 



(22) 



where n = (a^a)^ is the number of oscillator excitations corresponding to the instant magnetization. In the complete 
absence of the external field {ho ~ 0) Eq.([221) amounts to 

dn 1 

and the net magnetization relaxes from an initial value to the equilibrium one no = n(f2) = (e''^ — l) ^ with the 
spin-lattice relaxation rate T^^ = V^l- In general, by using the coefficients p^. we can recast Eg. ip^ as 



|J = (" - "0) - V2nm + hit) - hit), 



(23) 



where the time-dependent coefficients are 



_ V2hoiy'^ 
J^W - ^2 ,02 



{c.m<^' it) - n^" it)) , hit) 



(ri$'(t) + C™^"(i) 



As can be inferred from Eq. pD]) . /i,2(0 decreases exponentially with time and falls to zero as t > 
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The differential equation (|23|) can be solved numerically with the initial condition n{to) =0 {to < 0). The results 
for the simplest one-exponential case m = 1 are presented in Fig. 4 where the experimental data are plotted for 
comparison. We can see how the magnetization, following the h{t) variation, increases initially with time. After 
switching off the field, the magnetic moment of the system has a two-step evolution. The first, rapid stage ends when 
the system arrives at the critical state due to the balance of the magnetic driving force and the coercive force on the 
domain walls. The second, slow stage of the evolution is due to backward domain wall movement accompanied with 
damped oscillations around the moving center position. The damping indicates the effect of spin-lattice coupling. At 
rather low temperatures below 3D ordering, most of the system energy is lost to the bath because of dissipation. It 
is obviously seen from the inset of Fig. 4 that the damping is governed by the force fluctuation decay C in the bath. 
The slower the decay rate the more prominent domain wall oscillations. 

VI. CONCLUSIONS 

In conclusion, we summarize briefly the results presented in the paper. Measurements of the magnetization in 
pulsed fields for the molecule-based magnetic complex [Mn(hfac)2BN0H] composed from antiferromagnetically ordered 
ferrimagnetic chains show that the change in the magnetization in this compound below 3D ordering temperature 
is a slow dynamical process controlled by motion of magnetic domains with a single energy barrier height. As a 
critical field for the direct AF-FI transition is approached during application of a pulsed field, the domain wall motion 
along the separate field unfavorable chains starts to develop. After the field is switched off, the inverse transition 
from the field-induced FI state to the initial AF state in zero field provides very slow dynamics of the magnetization. 
The energy barrier hampering the reversal of the magnetization is originated from the antiferromagnetic interchain 
coupling. The latter thus plays an analogous role to that of magnetic anisotropy in a molecular cluster or a single 
chain magnet. The measurements of the ac-susceptibility in the region of the metamagnetic transition, where the AF 
and FI phase coexist, show that a leading contribution to the ac susceptibility results from the motion of domain 
walls separating the AF and FI phases. 

The domain-wall motion in both the ac and pulsed fields is accompanied by energy losses that causes a local 
heating of the samples. This is because the energy of the magnetic system transforms into the phonon energy, and, 
as a consequence, 3D magnetic ordering holds. Thus the system-bath coupling is a crucial in the description of the 
relaxation dynamics in [Mn(hfac)2BN0H]- 

On the basis of quantum dissipative theory in the standard second-order approximation for the system-bath Hamil- 
tonian, we derive Fokker-Plank equations for observables of the reference system. It is known that this well-justified 
approximation makes applicable the resultant FP equation in a large number of dissipative systems provided the 
system-bath coupling is not strong. The QDT-FP formalism has advantages of application convenience and straight- 
forwardness, as well as the ability to address both saturation processes caused by strong magnetic field pulses and 
processes far from a saturation by a small oscillating ac-field. 

The complex magnetic ac susceptibility is calculated within an approach, where the magnetization is controlled by 
the field-induced sideways motion of domain walls. The expressions for x'i'-^) ^-nd x"(<^) have the usual Debye form 
for small frequencies. ^From the maximum of the imaginary part of the ac susceptibility we evaluate the relaxation 
time that agrees qualitatively with that found from the remanent magnetization data. 

In the case of a small oscillating field the system remains near global equilibrium at all times, whereas a strong long- 
time driving field changing fast in comparison with relaxation to global equilibrium prohibits the normal evolution 
towards a Boltzmann distribution due to dynamical non-Markovian effects. In order to obtain a reliable understanding 
of the physics of the process we derive a FP equation in the framework QDT. The study of a time relaxation of a 
longitudinal magnetization shows that it experiences a two-step evolution after the field is switched off. The first 
rapid stage ends when the system arrives at the critical state where the magnetic driving force and the coercetive 
force acting jointly on the domain wall are balanced. The second slow stage of the evolution corresponds to backward 
domain wall movement together with damped oscillations around the moving domain-wall center. The damping is 
managed by a decay of force-force correlations of Langevin dissipative modes acting on the system from the bath. 

Acknowledgments 

We would like to thanks E.Z. Kuchinskii for discussions. This work was partly supported by the grant NREC-005 
of US CRDF (Civilian Research & Development Foundation). Two of us (V.E.S. and A.S.B.) thanks the Foundation 
Dynasty (Moscow) for the support. 



11 



VII. APPENDIX A. 



As a method for expressing the density operator in terms of c- number functions, the Wigner functions often lead 
to considerable simplification of the quantum equations of motion, as for example, by transforming operator master 
equations into more amenable Fokker-Plank differential equations. By the Wigner function one can express quantum- 
mechanical expectation values in form of averages over the complex plane (the classical phase space), the Wigner 
function playing the role of a c-number quasi-probability distribution 



f{z,z*;t) ^Tr p{t)Sw (a ~ z) 



where 



5w {a — z) 



(fx 



■ exp 



{ix* (a — z) + ix (a) — z*) } 



is the operator delta function with (Px — dxidx2, x = xi + ixi, z — zi + iz2. The Wigner function has the following 
property J f{z, z*- 1) d'^z = 1 and allows to easily evaluate expectations of symmetrically ordered products of the field 
operators, corresponding to Weyl's quantization procedure^^ 



where 



(z*rz"/(^,^*;t)rf'^ = ({(at)'"a"} 



and the symbol P denotes a permutation of the Bose operators. 
Using the last identity one obtains 



and 



y z* f{z, z*; t) d^ z = (a^)j , j z f[z, z*; t) d^ z = (a)^ 



z*zf{z,z*;t) d^z = (a^a)^ + 1/2. 



The Weyl symbol for any operator O is determined by 

(6) (z,z*)=TiTr 
V J w 

and the inversion formula is 



05w (a — z) 



0= I [0)^{z,z*)Swia- z)d^z. 



The Wigner functions of multiplication of two operators AB can be easily obtained from those of (^A^ and (^B^ 
using the following identities 



w 



V /w V /I 



z+l^,z*-l-^]{B) 



w \ 2du*' 2du 



w 



'w\ 2du* 



and Weyl symbols for Bose operators 

{a'')^{z,z*) = z\ {a)^{z,z*) 



41;) W.<"'°')i. 
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The relation between the Wigner function and the Weyl symbol of the density operator is the following 

f{z,z*;t) = 7r-'{ps{t))^{z,zn. 
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FIG. 1: Time dependence of the magnetization of [Mn(hfac)2BNOci] (a) and [Mn(hfac)2BNOH] (b) compounds at 1.5 K after 
application of the 5 T pulsed field. The logarithmic curve in the first case results from a relaxation mechanism that involves a 
distribution of energy barriers. The exponential decrease for the second case arises from a single-barrier activation mechanism. 



FIG. 2: In an external positive field the system consists of field- favorable (white) and field-unfavorable (shaded) chains. Interface 
kinks (domain walls) in the second-type chains can move only along the dotted arrows. 



FIG. 4: Model time dependence of the magnetization during and after application of the pulsed field found from the Fokker-Plank 
equation (dotted line). The experimental magnetization curve for [Mn(hfac)2BNOH] at T=1.61 K presented for comparison 
is shown as a solid line. The profile of the pulsed field used in the measurement is also plotted. The model parameters are 
v'^ = 0.3, = 0.95, T = 12, 7 = 0.004, Q, = 1.3, C,m ~ 0.7, ho = 1.3. Inset: time dependence of the longitudinal magnetization 
when 7 = 0.003, Cm = 0.4, = 0.35. The damped domain wall oscillations are resolved more clearly. 



FIG. 3: Temperature dependence of the relaxation time obtained from the remanent magnetization Mr at the value Mr = 1 A 
m^/kg after application of the 4 T pulsed field. The solid line shows Arrhenius behavior at lower temperatures. 



